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A THEORETICAL AND EXPERIMENTAL INVESTIGATION OF 
THE THREE-DIMENSIONAL VIBRATION CHARACTERISTICS OF A 
SCALED MODEL OF AN ASYMMETRICAL LAUNCH VEHICLE 

By Ellwood L. Peele, William M. Thompson, Jr., 
and Christine G. Pusey 
Langley Research Center 

SUMMARY 

The three-dimensional vibration characteristics of a dynamically scaled model of 
an asymmetrically clustered launch vehicle were investigated both analytically and experi- 
mentally. The matrix Holzer analysis formulated specifically for this model configura- 
tion is identified as being three-dimensional because it renders possible the analysis of 
motion which is coupled in three planes: lateral bending motion in two planes (pitch and 
yaw), torsional motion, and longitudinal motion. Previous two-dimensional analyses of 
similar configurations considered pitch-torsion uncoupled from yaw -longitudinal motion. 

The suitability of the analysis is evaluated by comparison of computed vibration 
characteristics with measured response of a near- replica scaled model to lateral excita- 
tion. On the basis of these comparisons, the procedure was found to be capable of pre- 
dicting the major vibration characteristics. Additional comparison is made between 
computed data based on a coupled and uncoupled mathematical model. These compari- 
sons indicate that for this configuration, the coupled and uncoupled analysis yielded essen- 
tially the same results. 

The study also demonstrates the value of using experimental data in conjunction 
with analytical results to obtain a keener insight into subtle structural behavior, and thus 
make a contribution to improved mathematical modeling. 

INTRODUCTION 

An understanding of the dynamic characteristics of launch vehicles is essential to 
insure reliable structural designs. Ideally, design data for complex launch vehicles 
should be based on the results of analytical computations; however, state-of-the-art 
methods are not always comprehensive or reliable enough to be the sole means of pre- 
dicting vehicle dynamic- response characteristics. Therefore, analytical methods for 
predicting the vibration response characteristics of structurally complex vehicles are 


usually verified by comparison with experimental data obtained from full-scale or 
dynamically scaled models. Several programs have been conducted with near-replica 
scaled models of complex launch vehicles at the Langley Research Center. These pro- 
grams include studies with models of launch- vehicle configurations consisting of: pro- 
pellant tanks symmetrically clustered around a multistage center tank (refs. 1 to 3), a 
tandem distribution of liquid propellant tanks in a single body (ref. 4), and propellant 
tanks asymmetrically clustered around a center body (ref. 5) . These studies have pro- 
vided experimental data which enabled verification of analytical techniques adapted to 
these configurations. 

The model configuration considered in this report represents a three-body asym- 
metrically oriented cluster. A matrix Holzer method has been formulated, in reference 6, 
for the calculation of the free vibration modes of multibodied launch vehicles both of the 
symmetrical and asymmetrical cluster configurations. Verification of the method was 
shown in reference 6 for the symmetric case through comparisons between calculated and 
measured vibration characteristics for the Saturn I model discussed in references 1 to 3. 
Similar comparisons for the asymmetric case, not given in reference 6, are given in this 
report. The purposes of the present report are (1) to demonstrate the applicability of 
the matrix Holzer fully coupled analysis to asymmetric clusters by comparisons between 
computed and experimentally determined modal data, (2) to compare theoretical computa- 
tions for both uncoupled and coupled analyses, and (3) to illustrate the use of experimental 
data both to improve the mathematical models employed in the analysis and to provide 
more accurate numerical input data. 

The experimental results were obtained during a ground vibration survey of a 
scaled dynamic model of the three-body clustered launch vehicle. For completeness, a 
brief review of the experimental phase of the test program together with a limited 
description of the model structure is presented in this report. Test results were obtained 
during lateral excitation of the model (pitch and yaw planes) when mounted in a suspension 
system designed to simulate free-flight boundary conditions. Data were obtained for dif- 
ferent propellant loading conditions simulating flight times which ranged from lift-off to 
solid- rocket-motor burnout and for three payload masses. 

SYMBOLS 


F longitudinal force 

f frequency (unmodified outrigger flexibility) 

f * frequency (modified outrigger flexibility) 
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acceleration due to gravity; structural damping 


area moment of inertia 
mass moment of inertia 
any dimensional length 

maximum model length, 257.6 in. (6.54 meters) 

yaw bending moment 

pitch bending moment 

element of generalized mass matrix 

mass 

torque 

time 

yaw shear force 
pitch shear force 

combined yaw-bending and shear deformation 
combined pitch-bending and shear deformation 
rectangular coordinates 
torsional angle of twist 
scale factor 

frequency parameter, co^/g 


longitudinal displacement 



cp slope of yaw-bending curve 

<p' slope of pitch-bending curve 

a) circular frequency, 2irf 

Subscripts: 

C beam component for center section 

c computed 

e experimental 

f full scale 

i * 1,2,3, . . . increment mass points of major beam components 
m model 

n number of lumped mass elements comprising beam component 

N beam component for nose section 

q discrete station along model span 

SRM1,SRM2 solid rocket motors 1 and 2, respectively 

T beam component for tail section 

u uncoupled 

I,n,m, . . . major internal boundaries 

Superscript: 

T transpose of a matrix 
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Matrices: 


O 

[] 


M.M 

[u<x)j 

[u(x)J 

<?}i 

< 2 } 





L°Jcluster 


column matrix 

square or rectangular 

transfer matrices defined in reference 6 

transfer matrix for entire model 

submatrix of |jJ(X)J 

state vector at ith station 

state vector of intermediate unknowns 

state vector with homogeneous boundary conditions imposed 

transfer matrix for ith elastic massless rod 

transfer matrix for lumped mass 

transfer matrix for mass and elastic rod combined 

submatrix (^crj^ compatible with {z} 

transfer matrix for SRM1, center core, SRM2 cluster 


ANALYSIS 


The analysis used to compute the model vibration characteristics is an application 
of the matrix Holzer or transfer- matrix method developed in reference 6. The matrix 
Holzer analysis, in review, is a matrix formulation of the tabular Holzer procedure which 
yields a solution, in stepwise manner, of the space-dependent part of a separable partial 
differential equation. The time-dependent solution is assumed sinusoidal. A set of fre- 
quencies are found, by trial and error, which give rise to inertial loadings that satisfy 
the external boundary conditions. This approach, also known as the transfer- matrix 
method, is discussed fully in reference 7. In reference 6, the procedure usually employed 
for single-beam configurations has been extended to provide the capability to accommodate 
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three-dimensional space-frame configurations. A complete derivation of the governing 
equations is presented in reference 6; therefore, only a brief outline of the analytical 
technique and its present application is given in this section. 

Model Configuration 

The subject of this investigation is a three-body configuration, shown in figures 1 
and 2, which consists of a center core containing simulated liquid propellants and two 
strap- on simulated solid-propellant rocket motors oriented diametrically opposite to one 
another at the aft end of the core. The plane of the core and solid- rocket motor center 
lines is defined as the yaw plane. The pitch plane is perpendicular to the yaw plane. 

The core consists of a cylindrically stiffened shell structure with fluid and mass 
elements distributed throughout its length. Basic components of the core are indicated 
in figures 1 and 2. The number of payloads and complexity of the stage III section 
requires additional comment. Two types of simulated payloads were employed: (1) the 
26 000 lbm (11 790 kg) and 45 000 lbm (20 410 kg) simulations which consist of a thick 
cylindrical shell bolted to the top of the core and (2) the 5000 lbm (2268 kg) payload 
(fig. 2) which consists of a rigid tubular shaft flexibly mounted on the top of the core. 

The stage III section is considered complex since mass eccentricities are introduced by 
the masses used to simulate propellants and tanks which are not oriented in a radially 
symmetric manner. 

The solid rocket motors (SRM) consist of relatively rigid steel cylinders which are 
attached by outriggers to the center core at the forward and aft ends as indicated in fig- 
ure 2. The forward outriggers are simple tie rods designed to transmit axial loads only. 
Therefore, pitch and yaw bending moments and longitudinal forces are not transmitted 
between the core and the solid rocket motors. The aft outriggers consist of trusses which 
are rigidly attached to the solid rocket motor aft sections. The core-truss attachment is 
made through two ball joints on each side. These joints are assumed to be capable of 
transmitting all forces and moments with the exception of yaw bending moments. A 
thrust- vector control tank, containing simulated fuel, is attached to each of the solid 
rocket motors as indicated in figures 1 and 2. 

Mathematical Model 

The space-frame mathematical model, shown schematically in figure 3, consists 
of the following beam components: nose section, center core, solid rocket motor 1, solid 
rocket motor 2, and the tail section interconnected by the forward and aft outriggers. 

Each component is free to bend in the pitch and yaw planes, to twist, and to undergo axial 
deformation subject to appropriate boundary conditions. The beams representing cylin- 
drical shell portions of the model are assumed to behave as massless rods having the 
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cross-sectional stiffness properties of the shell structure. The actual mass of the 
physical model structure is assumed to be concentrated at discrete points along each 
beam span and thus further subdivides the beam into elastic rod elements. The effect of 
mass eccentricity, which occurs when the center of gravity of a section is noncoincident 
with its centroid, is included. This condition exists in the model stage in because of the 
propellant-tank orientation and at the lower end of the solid rocket motors because of the 
thrust vector control tanks (TVC). (See fig. 2.) The thrust vector control tanks were 
assumed to be rigidly attached to the solid rocket motors and the mass was combined 
with the solid rocket motor mass. The addition of the mass eccentricity leads to coupling 
between pitch-torsion and yaw- longitudinal motion provided that the product- of- inertia 
terms between pitch and yaw are nonzero. The matrix Holzer analysis utilized herein is 
formulated to account for this fully coupled motion of the clustered configuration; this 
analysis is used to obtain the coupled response. However, the analysis may also be used 
to compute pitch-torsional motions separately from yaw -longitudinal motions and the 
results from this analysis are designated uncoupled responses. 

Liquid masses were assumed to be rigidly attached to the beam center line in the 
manner of structural masses. Fluid dynamics were assumed to be negligible. Rotatory 
and torsional moments of inertia of the liquids were assumed to be zero. 

Outline of Analytical Procedure 

The analytical development is basically a problem of obtaining a matrix relation- 
ship between the load and displacement quantities at some exterior boundary to those at 
another, that is, the transfer matrix. These variables, also described at interior points, 
consist of loads and displacements for pitch, yaw, torsion, and longitudinal motion. Since 
each mass point is assumed to be unrestrained to ground, it will have six degrees of free- 
dom and six corresponding loads. These quantities, given in vector form, are defined as 
the state vector (ref. 7) as follows: 

(z) = {V, M, <P,Y\ V', M', <*>', Y' ! F, £ i T, 0> T 

in which: 

V shear force 

M bending moment 

<p bending slope 

Y displacement 
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F longitudinal force 

£ longitudinal displacement 

T axial torque 

6 angle of twist 

Pitch and yaw directions are indicated by primed and unprimed symbols, 
respectively. The sign conventions established for these variables, when referred to 
any section of the model considered as a free body, are shown in figure 3(c). 

The exterior stations chosen for this analysis are the nose, station I, and the aft 
engine, station X. The transfer between these stations is accomplished by first deter- 
mining the state-vector transfer across elemental lumped-mass — weightless- rod line 
structures. The procedure, discussed in reference 6, is briefly as follows: the state 
vector is related to through the transfer matrix £Vj or 

( z )i + i = Hi( z >i <» 

where 

The matrices [ej i and jjtnj. are transfer matrices for the individual ith massless 
rod and concentrated mass, respectively. The matrices used in this paper are exten- 
sions of those presented in reference 6 and are developed to satisfy the need for a more 
detailed description of the stiffness distribution and the inclusion of open-ended branch 
beams. The matrices related to the beam elements are given in detail in the appendix. 
Once the transfer matrices, and Qn| ., for the rod- mass elements are obtained, 

the following relationships can be established for each beam component: 

®p,VI c ,VI 1> Vl2,Xj * [ a ](N,C,SRMl,SRM2,T)‘C Z }^I,V c ,V 1 ,V 2 ,IX) (2a) 

For example, the transfer across the nose section is given by 

{ z }n - [>]n©i < 2b > 

in which 

HN'HnHnCOn-lHn-l • • • [ f ]i[ m ]i ' ' ' [ £ ]l[ m ]l (2c) 

In equation (2c), n denotes the number of lumped mass elements comprising the nose 
section. 
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The transfer from station II across the entire three-beam cluster to station IX is 
then expressed by 


@ix - M 


SRM l+SRM2+C+outrigger 


“ r^cluster©!! 


(3) 


The derivation of the transfer matrix C°i) cluster is ver y involved and is not presented 
here. Details of the derivation are given in reference 6. In the derivation, the internal 
boundary conditions of continuity and equilibrium at the forward and aft outriggers are 
used to obtain solutions for certain intermediate unknowns at the solid rocket forward 
ends. The vector consisting of the intermediate unknowns is defined as 


Solutions for 
relation 


<*>v k -<v f V’,Tjv, V’, 


T j <t>, <*>', £;<*>, 4>\ |} y 

2 1 2 k 


(k= 1, 2) 



in terms of the state vector at station I are obtained from the 


®v k = M‘ 1 [ p ]( z >i 


(4) 


in which * and £pj are given in reference 6. Flexibility relationships for the 

forward outrigger and aft truss involved in the derivation of the and matrices 

are shown in the appendix. The intermediate unknowns are dependent upon the frequency 


parameter X = — assumed for the time-dependent motion and therefore must be deter- 
mined for each trial frequency. 


When the transfer relationship across the three-beam cluster, together with the 
nose and tail sections, is available, the entire model is spanned by multiplying all trans- 
fer matrices together to obtain 


{z} x - [UW] {z}j (5) 

in which [u(X)j represents the transfer matrix for the entire model. The remainder 
of the problem is the customary one of solving the set of homogeneous algebraic equations 

(6) 


where [lJ(X)^ is a submatrix of |u(X)j relating the external boundary conditions at X 

to those at I. The state vectors »x and < C z 3’j refect the homogeneous boundary 
conditions existing at these stations. The necessary and sufficient condition for a non- 
trivial solution for to exist is that the determinant jjJ(X)] = 0. This condition 

cannot be identically satisfied on digital computers but is approached numerically by 
trial- and- error methods to as many significant figures as is practical. In practice, the 
frequency parameter is varied over a wide range and a plot is made of the variation of 
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|U(X)J with A. Samples of these plots are shown in figure 4 for the matrix Holzer anal- 
ysis considering both a coupled and uncoupled mathematical model. Note that the plot is 
discontinuous for certain values of A. These values correspond to branch frequencies 
and frequencies for which [qJ in equation (4) is singular. In either case, the trans- 
fer matrices will become infinitely large. Approximate zeros of this function are 
obtained along with the components of the corresponding vector {z^. Deflection shapes 
are then constructed from the state vector components corresponding to each mass point. 
These components are obtained numerically by transferring the variables {V}j across 
the model span; for example, the state vector at station II (fig. 3(b)) is found from 

in which [o[)n * s a su t )ma t r i x of [ ct 3n compatible with The types of modal 

motions anticipated for the asymmetric cluster are illustrated in figure 5. 

Numerical Input 

The input data-used in the computations were largely obtained from the model manu- 
facturer. The various stiffness properties of the core and solid- rocket- shell structures 
are shown in figure 6. These parameters were obtained analytically and no experimental 
verification of these data exists. The mass distribution of the structure is given in 
table I which includes the center- of-gravity location and inertial properties. The masses 
representing liquids and distributed shell structure were further segmented into smaller 
lumped masses. The number of lumped masses varies for each weight condition; how- 
ever, an average of approximately 75 was employed for each configuration. Thus, sys- 
tems of 450 degrees of freedom were analyzed. 

Payload mounting fle xibility.- The differences in physical structure between the 
light and heavy simulated payloads require a different mathematical representation of 
each payload. The light payload (5000 lbm (2268 kg)) was interpreted as an elastically 
mounted branch beam; however, the complexity of the mounting prohibited an accurate 
analytical determination of the flexibility. The mounting flexibility was therefore based 
on values which are related to experimentally determined pitch and yaw frequencies of 
35.9 and 14.8 Hz, respectively. The heavier payloads (26 000 lbm (11 790 kg) and 
45 000 lbm (20 410 kg)) were each considered to be simple continuations of the main 
structure, that is, the core of the model having stiffness and mass properties as shown, 
respectively, in figure 6 and table I. 

Outrigger flexibilities .- Flexibility coefficients for the forward and aft outriggers 
were determined by the model manufacturer. These data, with indicated modifications, 
are presented in matrix form in table II and are defined by the equations shown in the 
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TABLE I.- MODEL COMPONENTS 


Items 

Mass 

Center of gravity 

Inertia 

X 

y 

z 

Roll 

Pitch 

Yaw 

Ibm 

kg | 

x/l 

in. j 

m ; 

in. | 

m 

AAlA 

m 

lb- in 2 

N/m 2 

lb- in 2 | 

N/m 2 j 

lb- in 2 

N/m 2 


Payloads 


Light 

35.7 

16.2 

-0.007 

-1.77 

-0.045 

0 

0 

12.00 

0.305 

531 

1.52 

2 310 

6.63 

2 310 

6.63 

Medium 

206.9 

93.85 

-.094 

-24.15 

-.613 

0 

0 

12.0 

.305 

27 901 

80.07 

42 890 

123.1 

42 890 

123.1 

Heavy j 

359.5 | 

163.1 | 

-.046 

-11.81 | 

-.300 

0 . 

0 

12.0 ! 

.305 

49 008 

140.6 

127 919 ( 

367.1 

127 919 

367.1 


Stage in components 


Guidance truss 

6.81 

3.09 

0.078 

20.1 

0.510 

0.1 

0.002 

4.7 

0.119 

77 

0.221 

20 

0.057 

54 

0.155 

Equipment truss 

5.44 

2.46 

.077 

19.8 

.503 

1.1 

.028 

19.7 

.500 

57 

.164 

11 

.032 

53 

.152 

Fuel tank (full) 

63.75 

28.92 

.124 

32.0 

.813 

-6.9 

-.175 

12.0 

.305 

302 

.867 

3 057 

8.77 

3 057 

8.77 

Oxidizer tank (full) 

124.35 

56.40 

.107 

27.5 

.699 

4.6 

.117 

12.0 

.305 

859 

2.47 

5 749 

16.50 

5 749 

16.50 

Control module skirt 

2.19 

.99 

.082 

21.0 

.533 

0 

0 

12.0 

.305 

316 

.907 

186 

.534 

186 

.534 

Propulsion module and helium sphere 

5.00 

' 2.27 

.109 

28.1 

.714 

0 

0 

12.0 

.305 

350 

1.00 

246 

.706 

109 

.313 

Truss and engines 

5.63 

2.55 

.146 

37.6 

.955 

0 

0 

12.0 

.305 

329 

.944 

352 

1.01 

389 

1.12 

Stage II components 

Forward oxidizer skirt 

7.36 

3.34 

0.178 

45.9 

1.17 

0 

0 

12.0 

0.305 

1 050 

3.01 

1 194 

3.43 

1 194 

3.43 

Forward oxidizer dome 

3.02 

1.37 

.214 

55.2 

1.40 

.4 

.010 

11.7 

.297 

267 

.766 

149 

.428 

149 

.428 

Aft oxidizer dome 

1.60 

.73 

.268 

69.1 

1.76 

0 

0 

11.6 

.295 

121 

.317 

72 

.207 

75 

.207 

Aft oxidizer skirt 

4.15 

1.88 

.264 

68.0 

1.73 

0 

0 

12.0 

.305 

590 

1.69 

370 

1.06 

370 

1.06 

Forward fuel skirt 

2.90 

1.32 

.310 

79.8 

2.03 

0 

0 

12.0 

.305 

412 

1.18 

226 

.649 

226 

.649 

Forward fuel dome 

3.16 

1.43 

.303 

78.0 

1.98 

-.4 

-.010 

12.3 

.312 

279 

.801 

156 

.448 

156 

.448 

Aft fuel dome 

2.86 

1.30 

.357 

91.9 

2.33 

0 

0 

11.7 

.297 

230 

.660 

129 

.370 

129 

.370 

Aft fuel skirt 

4.89 

2.22 

.357 

91.9 

2.33 

0 

0 

12.0 

.305 

691 

1.98 

454 

1.30 

454 

1.30 

Equipment truss 

3.94 

1.79 

.285 

73.4 

1.86 

.4 

.010 

16.8 

.427 

17 

.049 

4 

.011 

15 

.011 

Engine simulation 

11.44 

5.19 

.400 

103.1 

2.62 

0 

0 

12.0 

.305 

134 

.385 

208 

.597 

208 

.597 

Pressure test lines in oxidizer tank 

2.13 

.97 

.227 

58.4 

1.48 

-3.4 

-.086 

15.0 

.381 

83 

.238 

44 

.126 

39 

.126 

Pressure test lines in fuel tank 

2.13 

.97 

.316 

81.3 

2.07 

-3.4 

-.086 

15.0 

.381 

83 

.238 

44 

.126 

39 

.126 

Instrumentation in oxidizer tank 

1.00 

.45 

.227 

58.4 

1.48 

-3.4 

-.086 

15.0 

.381 

39 

.112 

21 

.060 

18 

.060 

Instrumentation in fuel tank 

.81 

.37 

.316 

81.3 

2.07 

-3.4 

-.086 

15.0 

CO 

co 

32 

.092 

17 

.049 

15 

.049 

Oxidizer fluid 

343.7 

155.90 

.242 

62.3 

1.58 

0 

0 

12.0 

.305 

0 

0 

18 688 

53.63 

18 688 

53.63 

Fuel fluid 

192.7 

87.4 

.330 

j 85.1 

2.16 

0 

0 

12.0 

.305 

0 

0 

8 608 

24.70 

8 608 

24.70 
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TABLE L- MODEL COMPONENTS - Concluded 


Items 

Mass 

1 

Center of gravity 


z 



Inertia 


iw 

J x 

i 


Etoll 

1 Pitch 

Yi 


Ibm 

) kg 

| x/l 

in. 

M 

in. 

M 

in. 

1 M 

lb-in^ 

| N-zn 2 

J lb-in 2 

| N-m2 

lb-in 2 

j N-m 2 





Stage I components 









Transportation section 

9.40 

4.26 

0.421 

108.4 

2.75 

0 

0 

12.0 

0.305 

1 369 

3.93 

1 310 

3.76 

990 

2.84 

Forward oxidizer skirt 

2.81 

1.27 

.492 

126.8 

3.22 

0 

0 

12.0 

.305 

399 

1.15 

213 

.611 

213 

.611 

Forward oxidizer dome 

3.08 

1.40 

.486 

125.1 

3.18 

0 

0 

12.0 

.305 

272 

.781 

152 

.436 

152 

.436 

Oxidizer tank barrel 

15.56 

7.06 

.594 

153.0 

3.8S 

.8 

.020 

12.8 

.325 

1 614 

4.63 

2 640 

7.58 

2 640 

7.58 

Aft oxidizer dome 

2.30 

1.04 

.709 

182.7 

4.64 

0 

0 

12.0 

.305 

166 

.476 

92 

.264 

92 

.264 

Aft oxidizer skirt 

2.47 

1.12 

.700 

180.4 

4.58 

0 

0 

12.0 

.305 

349 

1.00 

194 

.557 

194 

.557 

Forward fuel skirt 

4.06 

1.84 

.740 

190.7 

4.84 

0 

0 

12.0 

.305 

550 

1.68 

898 

2.58 

898 

2.58 

Forward fuel dome 

2.94 

1.33 

.746 

192.1 

4.88 

0 

0 

12.0 

.305 

260 

.746 

145 

.416 

145 

.416 

Fuel tank barrel 

18.90 

8.57 

.841 

216.7 

5.50 

.4 

.010 

12.4 

.315 

1 711 

4.91 

3 389 

9.73 

3 388 

9.73 

Aft fuel cone 

4.88 

2.21 

.945 

243.5 

6.18 

0 

0 

12.0 

.305 

281 

.806 

210 

.602 

210 

.602 

Aft fuel skirt 

7.31 

3.32 

.937 

241.4 

6.13 

0 

0 

11.9 

.302 

995 

2.86 

2 889 

8.29 

2 889 

8.29 

Engines 

37.73 

17.11 

1.000 

257.6 

6.54 

0 

0 

12.0 

.305 

1 099 

3.15 

3 060 

8.78 

3 060 

8.78 

Pressure test lines in oxidizer tank 

2.13 

.97 

.486 

125.1 

3.18 

-3.4 

-.086 

15.0 

.381 

83 

.238 

44 

.126 

39 

.112 

Pressure test lines in fuel tank 

2.13 

.97 

.746 

192.1 

4.88 

-3.4 

-.086 

15.0 

.381 

83 

.238 

44 

.126 

39 

.112 

Oxidizer fluid 

1321.1 

599.2 

.604 

155.5 

3.95 

0 

0 

12.0 

.305 

0 

0 

435 200 

1249 

435 200 

1249 

Fuel fluid 

689.5 

312.8 

.845 

217.8 

5.53 

0 

0 

12.0 

.305 

0 

0 

151 680 

435.3 

151 680 

435.3 




One solid rocket motor (empty) 








Nose section 

30.5 

13.8 

0.348 

89.56 

2.27 

25.3 

0.643 

12.0 

0.305 

2 027 

5.82 

2 595 

7.45 

2 595 

7.45 

Barrel 

270.4 

122.7 

.659 

169.81 

4.31 

25.3 

.643 

12.0 

.305 

38 601 

110.8 

477 556 

1370 

477 556 

1370 

Aft skirt 

53.19 

24.1 

.969 

249.60 

6.34 

25.3 

.643 

12.0 

.305 

7 048 

20.23 

4 498 

12.91 

4 498 

12.91 

Aft cone 

7.50 

3.4 

.958 

246.90 

6.27 

25.3 

.643 

12.0 

.305 

594 

1.70 

330 

.947 

330 

.947 

Engine 

58.63 

26.6 

1.000 

257.60 

6.54 

25.8 

.655 

12.0 

.305 

1 975 

5.67 

3 183 

9.13 

3 183 

9.13 

Forward outriggers 

1.16 

.53 

.390 

100.5 

2.55 

13.3 

.338 

12.0 

.305 

167 

.479 

167 

.479 

22 

.063 

Aft outriggers 

7.19 

3.26 

.952 

245.27 

6.23 

13.3 

.338 

12.0 

.305 

1 034 

2.97 

1 034 

2.97 

135 

.387 

Thrust vector control tank 

65.65 

29.8 

.721 

185.77 

4.72 

18.2 

.462 

25.9 

.658 

1 103 

3.17 

37 560 

107.8 

37 560 

107.8 

Barrel (full) 

305.69 

138.7 

.665 

171.31 

4.35 

25.3 

.643 

12.0 

.305 

43 638 

125.2 

539 883 

1549 

539 883 

1549 




One solid rocket motor (full) 








Nose section 

30.5 

13.8 

0.348 

89.56 

2.27 

25.3 

0.643 

12.0 

0.305 

2 027 

5.32 

2 595 

7.45 

2 595 

7.45 

Barrel 

305.69 

138.7 

.665 

171.31 

4.35 

25.3 

.643 

12.0 

.305 

43 638 

125.2 

539 883 

1549 

539 883 

1549 

Aft skirt 

53.19 

24.1 

.969 

249.60 

6.34 

25.3 

.643 

12.0 

.305 

7 048 

20.23 

4 498 

12.91 

4 498 

12.91 

Aft cone 

7.50 

3.40 

.958 

246.90 

6.27 

25.3 

.643 

12.0 

.305 

594 

1.70 

330 

.947 

330 

.947 

Engine 

58.63 

26.6 

1.000 

257.60 

6.54 

25.8 

.655 

12.0 

.315 

1 975 

5.67 

3 183 

9.13 

3 183 

9.13 

Forward outriggers 

1.16 

.53 

.390 

100.5 

2.55 

13.3 

.338 

12.0 

.305 

167 

.479 

167 

.479 

22 

.063 

Aft outriggers 

7.19 

3.26 

.952 

245.27 

6.23 

13.3 

.338 

12.0 

.305 

1 034 

2.97 

1 034 

2.97 

135 

.387 

Thrust vector control tank 

65.65 

29.8 

.721 

185.77 

4.72 

18.2 

.462 

25.9 


1 103 

3.17 

37 560 

107.8 

37 560 

107.8 





Ballast* ** per side 









Weight 1 

330.30 

149.8 

0.429 

110,52 

2.81 

25.3 

0.643 

12.0 

0.305 

24 624 

70.67 

12 577 

36.09 

12 577 

36.09 

*Weight 2 

329.10 

149.3 

.468 

120.52 

3.06 

25.3 

.643 

12.0 

.305 

24 335 

70.41 

12 532 

35.96 

12 532 

35.96 

Weight 3 

329.70 

149.5 

.538 

138.52 

3.52 

25.3 

.643 

12.0 

.305 

24 580 

70.54 

12 554 

36.03 

12 554 

36.03 

^Weight 4 

331.40 

150.3 

.577 

148.52 

3.77 

25.3 

.643 

12.0 

.305 

24 706 

70.90 

12 614 

36.21 

12 614 

36.21 

Weight 5 

330.40 

149.9 

.650 

167.52 

4.26 

25.3 

.643 

12.0 

.305 

24 632 

70.69 

12 581 

36.11 

12 581 

36.11 

^Weight 6 

331.50 

150.4 

.685 

176.55 

4.48 

25.3 

.643 

12.0 

.305 

24 714 

70.92 

12 623 

36.23 

12 623 

36.23 

Weight 7 

330.80 

150.0 

.763 

196.55 

4.99 

25.3 

.643 

12.0 

.305 

24 662 

70.78 

12 596 

36.15 

12 596 

36.15 

^Weight 8 

329.80 

149.6 

.802 

206.55 

5.25 

25.3 

.643 

12.0 

.305 

24 587 

70.56 

12 558 

36.04 

12 558 

36.04 

Weight 9 

329.60 

149.5 

.879 

226.55 

5.75 

25.3 

.643 

12.0 

.305 

24 572 

70.52 

12 551 

36.02 

12 551 

36.02 

^Weight 10 

330.40 

149.9 

.918 

236.55 

6.01 

25.3 

.643 

12.0 

.305 

24 632 

70.69 

12 581 

36.11 

12 581 

36.11 

**Thrust vector control propellant (1/2 full) 

62.0 

28.1 

.894 

230.41 

5.85 

18.2 

.462 

25-9 

.658 


— 

2 855 

8.19 

2 855 

8.19 

**Thrust vector control propellant (full) 

123.90 

56.2 

.851 

219.20 

5.57 

18.2 

.462 

25.9 

.658 


— 

21 249 

60.98 

21 249 

60.98 


*Omit to obtain data for one solid rocket motor (1/2 full). 

**Omit to obtain data for one solid rocket motor (full). 
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TABLE n.- OUTRIGGER FLEXIBILITY COEFFICIENTS 


[The coefficient matrices are defined by equations appearing in the appendixj 


1.6187 

0.0777 

0.5674 

0.0292~ 


2.0663 

0.0444 

0.3240 

0.0167 

0.0777 

0.003375 

0.0292 

0.00155 

x 10" 5 in./lb = 

0.0444 

0.001927 

0.0167 

0.000885 

0.5674 

0.0292 

3.6187 

0.0777 

0.00324 

0.0167 

2.0663 

0.0444 

0.0292 

0.00155 

0.0777 

0.0034 


0.0167 

0.000885 

0.0444 

0.001927 


[B> 


26.15 

6.9 


6.9 

26.151 


X 1<T 6 in./lb 


14.93 

3.9 


3.9 

14.93 


X 10" 8 


m/N 


H= 


2.8172 

0.0428 

0.01695 

0.6537 

0.02115 

0 

1.6086 

0.02444 

0.009679 

0.3733 

0.01208 

0 

0.0428 

0.00717** 

0.00015 

0.02115 

0.00478** 

0 

0.0244 

0.00409** 

0.00086 

0.01208 

0.00273** 

0 

0.01695 

0.00015 

0.040* 

0 

0 

0 

x 10"' 5 in./lb = 

0.009679 

0.000086 

0.023* 

0 

0 

0 

0.6537 

0.002115 

0 

2.8172 

0.0428 

0.01695 


0.3733 

0.01208 

0 

1.6086 

0.0244 

0.009679 

0.02115 

0.00478** 

0 

0.0428 

0.00717** 

0.00015 


0.01208 

0.00273** 

0 

0.0244 

0.00409** 

0.000086 

0 

0 

0 

0.01695 

0.00015 

0.040* 


0 

0 

0 

0.009679 

0.000086 

0.023* 


M- 


0 

1 

0 

0 

0 

0 

2.575 X 10" 6 (in./lb) 

6.75 x lO -6 (in./lb) 


0 

1 

0 

0 

0 

0 

1.470 x 10' 8 (m/N) 

1 2 

£ 

00 

1 

o 

X 

CO 

■'*. 

o 

l - a 

0 

0 

0 

0 

1 

0 

0 

= 

l - a 

0 

0 

0 

0 

1 

0 

0 

0 

1 

0 

0 

0 

0 

0.75 x 10- 6 (in./lb) 

2.575 x lO’ 6 (in./lb) 


0 

1 

0 

0 

0 

0 

0.43 x 10“ 8 (m/N) 

1.470 x 10- 8 (m/N) 

-(1 - a) 

0 

0 

0 

0 

1 

0 

0 


-a - a) 

0 

0 

0 

0 

1 

0 

0 


6.9 

-2.55 

0 

0 


3.9 

-1.46 

0 

0 

-2.55 

2.10 

0 

0 

x 10' 6 in./lb = 

-1.46 

1.20 

0 

0 

0 

0 

6.9 

2.55 

0 

0 

3.9 

-1.46 

0 

0 

-2.55 

2.10 


0 

0 

-1.46 

1.20 
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^Revised to reflect assumption of negligible pitch-bending moment restraint. 
**Re computed to be consistent with yaw-bending coefficients. 
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appendix. The coefficients reflect the assumptions that the forward outriggers are 
incapable of transmitting pitch and yaw bending moments or axial force. The aft out- 
rigger coefficients reflect the assumptions that all forces and moments except yaw 
bending can be transmitted. Outrigger coefficients computed under these assumptions 
are identified as unmodified. During the present investigation, comparisons between 
computed and experimental vibration test results led to a reappraisal of the assumption 
that the aft outriggers could transmit pitch bending moment between the solid rocket 
motors and the core. The reasons for the reappraisal and the modifications in the coef- 
ficients are discussed fully in the section "Results and Discussion." Modal data uti- 
lizing both modified and unmodified outrigger coefficients were computed. 

EXPERIMENTAL TEST APPARATUS AND PROCEDURE 
Model- Scaling Considerations 

Only a general description of the overall structure is presented in this report since 
a more detailed description of the model center core is available in reference 4. The 
materials used in the structure, primarily aluminum alloy, are the same as those used 
for the full-scale vehicle of which the model is an approximate structural replica. In 
some noncritical areas, different alloys of aluminum were substituted to facilitate fabri- 
cation. Whenever practical, geometric scaling was used in the design of the model to 
insure faithful reproduction of the main load paths and thus lead to true dynamic charac- 
teristics. It should be emphasized, however, that the structure is not a true replica 
model in that some components were represented by simple structural designs which are 
only dynamically similar to the prototype. An example of this type of dynamic scaling 
in basic structural members is the substitution of simple shapes for shell stiffeners for 
which the cross-sectional area and area moment of inertia were scaled. Other examples 
include lumped mass representations of equipment items, engines, and, in some cases, 
propellants. The various structural parameters of the model were scaled by the fol- 
lowing relationships in which r] is the scale factor: 

Length: 

^m = 

Area moment of inertia: 

im = 

Mass: 

m m = *? 3m f 

Frequency: 

w m * \ w f 

Mass moment of inertia: 

I m = ?? 5 I f 
14 


Model Description 


The center core of the model is a sectional structure approximately 23 feet (7 m) 
in length consisting of a payload, stage HI section, stage II, and stage I as major com- 
ponents. (See fig. 1.) 

The stage III section indicated in figure 1 is a relatively compact and massive com- 
ponent of the core in which such elements as engines and liquid propellant tanks were 
simulated by cylindrical steel weights. These masses were mounted on trusses in the 
stage in section so that the rigid-body dynamic properties simulated this part of the 
prototype. 

The stage II section consists of the fuel tank, oxidizer tank, tank skirts, and the 
propulsion system. The propellant- tank skirts are basically thin-skinned (0.009 inch 
(0.299 mm)) shells reinforced by longitudinal and circumferential stiffeners riveted to 
the skin. The oxidizer tank and the fuel tank were liquid loaded and pressurized during 
the vibration tests. The full-scale liquid propellants were simulated in the model by 
liquids having the same density as the actual propellants. The various densities are 
shown in table III. The single-engine propulsion system for stage II was simulated by a 
steel cylinder. 


TABLE HI.- SIMULATED LIQUID PROPELLANTS 


Propellant 

Mixture (by weight) 

Specific gravity 

Oxidizer 

50.2 Freon; 49.8 methylene chloride 

1.42 

Fuel 

57.0 alcohol; 43.0 water 

.93 


The transportation section (which acts as the interstage between stages I and II), the 
two large liquid propellant tanks, and the twin- engine propulsion module constitute the 
first stage of the launch vehicle. The transportation section is a thin-skinned (0.009 inch 
(0.229 mm)) longitudinally and circumferentially stiffened structure similar in construc- 
tion to the stage II tank skirts. The oxidizer and fuel tanks are also capable of pressur- 
ization and liquid loading with the same fluids (identified in table III). The propellant 
tanks in this stage are internally stiffened by longerons welded to the tank wall and cir- 
cumferential ring stiffeners welded to the longerons. 

The two identical solid rocket motors and the attached thrust vector control (TVC) 
fuel tanks compose the "zero stage" of the three-body configuration. Each solid motor 
case consists of two sectional steel barrels, a nose cone, and an aft skirt. The top and 
bottom of each solid rocket motor is attached to the center core by forward and aft out- 
riggers, respectively, as indicated in figure 2. The zero-stage solid fuel was simulated 
by lead disks (each disk had a nominal mass of 330 lbm (149.7 kg)) evenly distributed 
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throughout the length of the barrel section and attached in order not to affect stiffness. 
Three loading conditions of the solids - empty, one-half full, and full - were simulated 
by installing 0, 5, or 10 lead disks, respectively. Each of the associated thrust vector 
control (T VC) tanks consists of two attached cylindrical sections (see fig. 1) which were 
designed to be fluid loaded but not pressurized. 

A detailed listing of the mass and center-of-gravity location (normalized at sta- 
tion 257.6 in. (6.54 m)) of items comprising the major components of the model described 
is given in table I. The mass data for the model shown in table I were calculated by use 
of the dimensions and material density of the structural elements of the model. The mass 
of each payload is a measured value obtained by weighing the payload configuration pro- 
vided for the tests. The stiffness characteristics of the model were provided by the manu- 
facturer and are presented in figure 6. 

The three payloads representing masses of 5000 lbm (2268 kg), 26 000 lbm 
(11 790 kg), and 45 000 lbm (20 410 kg) full scale were not models of any known prototype 
design but provided a means to study the effect of payload mass on the vibration response 
of the model structure. 


Suspension System 

During the vibration tests, the model configuration was suspended in a cable and 
spring system which offered negligible restraint against oscillatory motion in the lateral 
planes. Rigid-body frequencies of the model-cable system were at least one decade 
below the lowest structural modes. The suspension, shown schematically in figure 7, 
consisted of cable supports on both sides of the model. On each side of the structure one 
end of the cable was attached directly to the solid rocket motor casing near the level of 
the forward shear ties. The other end passed through an overhead pulley and terminated 
at a leaf spring located near the base of the model. The model elevation was adjustable 
with the use of hand winches placed at the leaf spring end of the cable to compensate for 
elongation of the support cables caused by different fuel loading conditions. 

Shaker System 

Excitation of the model was accomplished by use of a set of eight matched electro- 
magnetic shakers each rated at 50 vector force pounds (222.4 N). The shakers were 
suspended from overhead by cables to eliminate support resonances and to permit nearly 
unrestricted choice of the location of shaker attachment points on the model. All shakers 
were attached to the model by a vacuum-type attachment which also permitted versatility 
in choosing shaker attachment locations. A diagram showing the total number of shaker 
locations for the model in the pitch and yaw planes is shown in figure 8. Some of the sus- 
pended shakers used to excite the model in the yaw plane may be seen in figure 2. Each 
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shaker was powered by an individual amplifier, but frequency of excitation on all shakers 
was simultaneously controlled by a common oscillator. Operation of the shaker system 
was centralized at a control panel which allowed the use of one or more shakers at the 
same time. The control panel was equipped with four dual-beam oscilloscopes on which 
the armature voltage and current of any shaker could be monitored. A switching arrange- 
ment was provided so that on any scope the phase of the voltage or current of a shaker 
armature could be observed relative to the phase of the velocity or current of the arma- 
ture of any other shaker. This arrangement enabled a check of the phase relationship 
between the different shakers used in exciting a particular mode. The control panel also 
contained phase switches which could change the phase of each shaker by 180°. Finally, 
a master switch was installed which allowed a simultaneous termination of power to all 
shakers in operation. This switch was used to remove the excitation force to obtain 
records of the decay of the model oscillation during the transient phase of the response. 

Instrumentation 

The mode shapes and frequencies at resonances of the model were determined with 
the use of two movable probes. The probes were crystal-type accelerometers coupled 
with a cathode follower and amplifier as signal conditioning equipment. In addition, a 
crystal-type force gage was linked in the arm of the shaker attachment at x/l = 0.934 
in order that the shaker excitation force could be measured. A schematic diagram of 
this instrumentation system is shown in figure 9. The accelerometers and force gage 
were calibrated through their respective signal conditioning systems and were found to 
have a flat response from 6 Hz to above 100 Hz. 

Test Procedure 

Oscillation of the model in the pitch and yaw directions was induced to obtain 
resonant frequencies and associated mode patterns for the test condition shown in table IV. 

TABLE IV.- SUMMARY OF MODEL TEST CONDITIONS 


Nominal payload 

Flight time, 
sec 

Excitation direction 

lbm 

kg 

5 000 

2 268 

0, 53, 105 

Pitch, yaw 

26 000 

11 790 

0, 105 

Pitch, yaw 

45 000 

20 410 

0, 53, 105 

Pitch, yaw 


In tuning the modes for different fuel loading conditions of the model, several combina- 
tions of shakers were required; however, the following description reviews the general 
procedure followed in exciting modes for all loading conditions. The shaker in the 
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location estimated to provide the maximum model excitation was used to excite the mode. 
The Lissajous figure for the variation of shaker armature current with voltage (variation 
of velocity with force) was displayed on the monitoring oscilloscope and the frequency 
adjusted until the current and voltage were in phase. Traces of the decay of the model 
deflection response were recorded in order that the decay frequency could be compared 
with the forcing frequency to ascertain accurate tuning of the mode. Other needed 
shakers were sequentially energized and shaker amplitude, polarity, and frequency were 
adjusted to insure that the shaker force was in phase with velocity at all points of excita- 
tion on the model. In addition, for the shaker located at x/l * 0.934 where the force 
gage was linked in the shaker arm, the variation of shaker force with velocity was also 
monitored. 

Mode shapes were measured with the use of a fixed and a movable probe by noting 
the amplitude and phase at various points on the model relative to a chosen antinode loca- 
tion for the particular mode. The amplified signals from these accelerometers were 
filtered to eliminate signal noise and were fed into an oscilloscope displaying a Lissajous 
figure normalized at a 45° slope at the location of the fixed accelerometer. A transparent 
circular plastic face' calibrated in tangents and attached to the front of the oscilloscope 
screen permitted direct reading of the amplitude and phase. Therefore, by moving the 
accelerometer to different locations on the model, the amplitude and phase of the response 
at any location on the model relative to the location of the fixed accelerometer was read 
directly from the scope face. 


RESULTS AND DISCUSSION 

The measured and computed vibration characteristics are presented for representa- 
tive payloads and flight conditions. Experimental results typical of the launch-vehicle 
response are compared with computed results and specific areas of disagreement are 
discussed in detail. To avoid unnecessary repetition of data, experimental and computed 
mode shapes are presented in the same figures. All mode shape data are shown in fig- 
ures 10 to 17 and associated frequencies in table V; however, the vibration data are 
discussed in the appropriate order under the following headings. 

Experimental Modes 

The experimental mode shapes and associated resonant frequencies were obtained 
for the test conditions shown in table IV for the tank pressures and propellant masses 
given in table VI. The experimental mode data presented in this paper are limited to 
typical examples of lateral response which represent the first three lateral modes for 
all propellant loading conditions and payload masses. These typical examples of mea- 
sured mode shapes are shown for the intermediate weight 26 000 lbm (11 790 kg) payload. 
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TABLE V.- COMPARISON OF MEASURED AND COMPUTED FREQUENCIES 








Frequency, Hz, for flight times of - 










0 sec for — 





53 sec for - 




105 sec for 

- 


Type of frequency 


Pitch 



Yaw 



Pitch 



Yaw 



Pitch 


Yaw 


Mode 

Mode 

Mode 

Mode 

Mode 

Mode 

Mode 

Mode 

Mode 

Mode 

Mode 

Mode 

Mode 

Mode 

Mode 

Mode 

Mode 

Mode 


1 

2 

3 

1 

2 

3 

1 

2 

3 

1 

2 

3 

1 

2 

3 

1 

2 

3 





Payload of 5000 lbm 

(2268 kg) 










Experimental modes, f e 

8.7 

13.0 

31.1 

8.4 

13.2 

21.7 

9.1 

14.2 

26.7 

9.3 

14.5 

20.1 

10.3 

19.8 

30.8 

10.6 

18.7 

24.2 

Coupled modes utilizing unmodified 

8.9 

16.2 

39.6 

7.4 

13.3 

15.9 







10.6 

23.3 

33.9 

9.2 

14.9 

22.4 

outrigger coefficients, f c 



















Coupled modes utilizing modified 

8.2 

13.6 

37.0 

7.4 

13.2 

15.9 

8.6 

14.7 

28.5 

8.7 

14.4 

21.6 

9.6 

20.6 

32.8 

9.2 

14.9 

22.5 

outrigger coefficients, f c * 



















Uncoupled modes utilizing modified 







8.6 

14.7 

28.5 

8.2 

14.4 

21.6 







outrigger coefficients, f u 























Payload of 

26 000 lbm 

(11 790 kg) 









Experimental modes, f e 

6.0 

12.7 

20.3 

5.8 

12.0 

18.9 







7.5 

16.2 

25.8 

7.5 

15.4 

24.5 

Coupled modes utilizing unmodified 

6.2 

15.6 

20.7 

5.5 

12.2 

21.4 







7.7 

18.3 

31.6 

7.1 

15.9 

30.1 

outrigger coefficients, f c 



















Coupled modes utilizing modified 

5.9 

12.6 

20.6 

5.5 

12.2 

21.4 







7.3 

15.9 

29.3 

7.1 

15.9 

30.1 

outrigger coefficients, f c 



















Uncoupled modes utilizing modified 

5.9 

12.6 

20.7 

5.5 

12.8 

21.4 







7.2 

15.9 

29.3 

7.1 

15.9 

30.1 

outrigger coefficients, f u 























Payload of 

45 000 lbm 

(20 410 kg) 









Experimental modes, f e 

5.6 

12.7 

20.3 

5.4 

12.0 

18.4 

6.1 

13.6 

20.6 

6.2 

13.3 

19.2 

7.3 

16.2 

25.9 

7.3 

15.5 

24.2 

Coupled modes utilizing unmodified 

5.8 

15.6 

20.5 

5.1 

12.2 

21.2 







7.3 

18.2 

31.2 

6.8 

15.8 

29.7 

outrigger coefficients, fc 



















Coupled modes utilizing modified 

5.5 

12.5 

20.5 

5.1 

12.2 

21.2 

6.0 

13.3 

21.3 

5.6 

13.2 

20.7 

7.0 

15.8 

27.1 

6.8 

15.8 

29.7 

outrigger coefficient, f c * 



















Uncoupled modes utilizing modified 

— 

— 

— 

___ 

— 

— 

6.0 

13.3 

21.3 

5.6 

13.2 

20.7 







outrigger coefficients, f u 




















TABLE VI.- TANK PRESSURE AND PROPELLANT LOADINGS 
(a) Tank pressure 





Tank pressure* for - 




Flight 


Oxidizer 



Fuel 


time, 

sec 

Stage I 

Stage n 

Stage I 

Stage II 


lb/in2 

N/ cm2 

lb/cm2 

N/cm 2 

lb/cm 2 

N/cm 2 

lb/in 2 

N/cm 2 

105 

7.0 

4.8 

9.0 

6.2 

8.2 

5.7 

9.4 

6.5 

53 

7.0 

4.8 

9.0 

6.2 

8.2 

5.7 

9.4 

6.5 

0 

7.0 

4.8 

9.0 

6.2 

8.2 

5.7 

9.4 

6.5 


(b) Propellant loadings 






Solid 

and liquid propellant loading for - 




Flight 

time, 

sec 


Stage n 



Stage I 


Stage 0* 

Thrust 

vector 

control** 

Oxidizer 

Fuel 

Oxidizer 

Fuel 

Solid** 

Tank** 


lbm 

kg 

lbm 

kg 

lbm 

kg 

lbm 

kg 

lbm 

kg 

lbm 

kg 

105 

343.7 

155.9 

192.7 

87.4 

1321.1 

599.2 

689.5 

312.8 

0 

0 

0 

0 

53 

343.7 

155.9 

192.7 

87.4 

1321.1 

599.2 

689.5 

312.8 

1650.0 

748.4 

62.0 

28.1 

0 

343.7 

155.9 

192.7 

87.4 

1321.1 

599.2 

689.5 

312.8 

3300.0 

1496.9 

129.3 

58.6 


*Only stage I and stage II tanks require pressurization. 
**Above loadings are for individual tanks. 
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The pitch and yaw responses for lift-off (t = 0) are shown in figures 12 and 13, respec- 
tively; and the asymmetric solid rocket motor pitch response is shown in figure 14. 
(Modes are considered symmetric if the motions of the solid rocket motors are in phase 
as indicated in fig. 5.) Corresponding lateral mode response at solid rocket motor burn- 
out (t = 105 sec) is shown in figures 15 to 17. In general, the motions for symmetric 
modes of the two solid rocket motors were coincident as indicated by the experimental 
deflection patterns in figures 12, 13, 15, and 16. At numerous stations on the model core, 
the displacements (figs. 12, 13, 15, and 16) were measured on the front (circular symbol) 
and back (square symbol) of the structure to indicate localized response and thus avoid 
giving a distorted view of the beam-type lateral motion. 

Theoretically, lightly damped structures at resonance have mode shapes which are 
orthogonal. ^Structural damping ratios observed during decay tests in the model were 
very low in the order of 1 percent critical damping; that is, g = 2 ^ = 0.02.^ The accu- 
racy of the symmetric experimental mode shapes is, therefore, indicated by computation 
of the generalized masses by utilizing all modes including rigid body translation and 
rotation. The generalized masses obtained by using the measured mode shapes pre- 
sented in this report were computed from the mass data given in table I and are shown in 
tables VII and VIII as normalized generalized mass matrices. The elements of the gen- 
eralized mass matrix M^j were calculated by use of the equation: 


N 

M ij = I yi( x q)yj( x q) m q = M ii 

N 

M ij = £ yi(x q ) yj (x q )m q -0 
q=l 

where 

m q mass located at station x q 

yi(x q ) ith modal displacement function evaluated at x q 

y j( x q ) jth modal displacement function evaluated at x q 


(i = j) 

> ( 8 ) 

(i* j) 
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TABLE VC.- ORTHOGONALITY CHECK OF EXPERIMENTAL MODE SHAPES 
[t = 0, 26 000 lbm (11 790 kg) payload] 

(a) Pitch (b) Yaw 


Elements of generalized mass matrix 


Elements of generalized mass matrix 

Mode 

Frequency, 

Mu 


Mode 

Frequency, 

Mii 

Hz 

lbm 

kg 


Hz 

lbm 

kg 



Rigid body translation 

0 

11 032.0 

5004.0 


Rigid body translation 

0 

11 032.0 

5004.0 

Rigid body rotation 

0 

962.2 

436.4 


Rigid body rotation 

0 

962.2 

436.4 

Mode 1 

6.0 

368.9 

167.3 


Mode 1 

5.8 

277.8 

126.0 

Mode 2 

12.7 

1 941.1 

880.5 


Mode 2 

12.0 

860.0 

390.1 

Mode 3 

20.3 

543.7 

246.6 


Mode 3 

18.9 

775.9 

351.9 


Normalized generalized mass matrix Normalized generalized mass matrix 


1.000 

0.002 

-0.028 

0.044 

-0.012 


1.000 

0.002 

-0.021 

0.027 

0.006 

0.002 

1.000 

0.019 

0.032 

-0.032 


0.002 

1.000 

0.068 

-0.001 

-0.008 

-0.028 

0.019 

1.000 

-0.102 

-0.034 


-0.021 

0.068 

1.000 

0.134 

-0.118 

0.044 

0.032 

-0.102 

1.000 

0.073 


0.027 

-0.001 

0.134 

1.000 

0.184 

-0.012 

-0.032 

-0.034 

0.073 

1.000 


0.006 

-0.008 

-0.118 

0.184 

1.000 


TABLE VUE.- ORTHOGONALITY CHECK OF EXPERIMENTAL MODE SHAPES 
[t = 105 seconds, 26 000 lbm (11 790 kg) payload] 


(a) Pitch 


Elements of generalized mass matrix 


Mode 

Frequency, 

M 

ii 

Hz 

lbm 

kg 

Rigid body translation 

0 

4123.4 

1870.3 

Rigid body rotation 

0 

700.4 

317.7 

Mode 1 

7.5 

383.9 

174.0 

Mode 2 

16.2 

1134.1 

514.4 

Mode 3 

25.8 

937.6 

425.3 


Normalized generalized mass matrix 


1.000 

-0.001 

-0.083 

0.012 

0.114 

-0.001 

1.000 

0.064 

-0.022 

0.068 

-0.083 

0.064 

1.000 

-0.070 

0.052 

0.012 

-0.022 

-0.070 

1.000 

-0.099 

0.114 

0.068 

0.052 

-0.099 

1.000 


(b) Yaw 


Elements of generalized mass matrix 


Mode 

Frequency, 

Mii 

Hz 

lbm 

kg 

Rigid body translation 

0 

4123.4 

1870.3 

Rigid body rotation 

0 

700.4 

317.7 

Mode 1 

7.5 

325.0 

147.4 

Mode 2 

15.4 

1152.9 

522.9 

Mode 3 

24.5 

633.0 

287.1 


Normalized generalized mass matrix 


1.000 

-0.001 

-0.047 

-0.046 

0.124 

-0.001 

1.000 

0.038 

-0.073 

-0.030 

-0.047 

0.038 

1.000 

0.021 

-0.038 

-0.064 

-0.073 

0.021 

1.000 

-0.028 

0.124 

-0.030 

-0.038 

-0.028 

1.000 
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The numbers tabulated under the heading "Elements of the generalized mass matrix" 
in tables VII and Vm are the Mjj terms computed from equations (8) and the model mass 
data and the measured mode shapes. The element of the normalized generalized mass 
matrix Mjj is defined by 

- M ij 
Mi] = = J 


Analytical Modes and Comparison With Experiment 

At the beginning of this section a selected number of computed lateral modes which 
incorporate unmodified aft outrigger coefficients are presented. These computed data 
are compared with experimental results to indicate the need for modification of the aft 
outrigger joint fixity. If the modified aft outrigger coefficients in the coupled mathemati- 
cal model are utilized, lateral modal data can then be presented for all propellant loading 
conditions and payload configurations. Longitudinal and torsional modes also obtained by 
using the modified aft outrigger coefficients in the coupled analysis are given for several 
representative cases. Additional computations (designated the uncoupled analysis) are 
discussed which show the effect of omitting inertial coupling due to mass eccentricity in 
the mathematical model. 


Effect of ou trigger flexibility on lateral modes.- Initially, analytical modes and 
frequencies were computed by using the unmodified flexibility coefficients of the aft out- 
riggers identified previously. Typical examples of the lateral mode shapes, thus 
obtained, are compared with measured data in figures 10 and 11. Computed frequencies 
f c are compared with measured values f e for the two extreme propellant loading con- 
ditions in table V. In these comparisons of computed and measured mode shapes and 
frequencies, disagreement is significant in the second mode. The source of the discrep- 
ancy was traced to causes exhibited by two interesting features in the measured and com- 
puted mode shapes: The first is the difference in the slope of the core displacement pat- 
tern at the aft end of the structure and the second is the difference in the displacements 
of the solid rocket motors. (See figs. 10 and 11.) In these initially computed curves, the 
slope at the aft end of the core ( x/l = 0.942) is in agreement with the slope of the solid 
rocket motion which reflects the assumed fixed restraint at the a ft outrigger attachment 
against rotation of the solid motors in the pitch plane ^that is, coefficients of Myjj 

and SRM2) & iven * n the appendix. The measured second mode shapes (figs. 10(b) 

and 11(b)), however, show the slopes of the core and solid motors to be nearly perpendic- 
ular and indicate that the fixed restraint used in the computation may be a misrepresenta- 
tion of the actual behavior of the model at this point. (Note that the end experimental 
point in fig. 11(c) denotes motion of the stage I engine and does not indicate a reversal of 
slope of the aft end of the model core.) 
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The vibration data were, therefore, recomputed by using modified coefficients 
which assume a negligible restraint, approaching a pinned condition, against rotation of 
the solid motors in the pitch plane (that is, My m gRMl = M vxn SRM2 ~ • These data 

constitute the major results of the analytical effort and are presented and discussed in 
this section. 

Lateral modes, coupled analysis .- Mode shapes were computed for all propellant 
loading conditions listed in table IV and for three payload masses; however, only a repre- 
sentative set having dominant motion in the pitch and yaw planes is included since a com- 
parison of these results indicates, in general, the degree of agreement attained. This 
set of mode shapes is shown in figures 12 to 17 and were obtained for the 26 000 lbm 
(11 790 kg) payload at flight times t = 0 and t = 105 seconds. The computed displace- 
ments are compared with measured results on the basis of relative amplitude and loca- 
tion of nodes and are generally in satisfactory agreement with the measured results. As 
evident in figures 12 to 17, the results compare more favorably for the fundamental mode 
but compare progressively less favorably for higher modes. Poor agreement in com- 
puted and measured mode shapes occurred in only three of the 48 computations: the 
pitch and yaw third mode of the 5000 lbm (2268 kg) payload, t = 0 flight time, and the 
third yaw mode of the 45 000 lbm (20 410 kg) payload, t = 53 seconds flight time. 

These cases, although not shown, possessed qualitative similarity in computed and mea- 
sured mode shapes but exhibited poor agreement in light of the amplitude-node location 
criteria used as a basis for comparison. 

The computed natural frequencies associated with the mode shapes discussed are 
compared with measured resonant frequencies and are shown for all flight conditions and 
payload weights in table V. Computed frequencies utilizing modified outrigger coeffi- 
cients are designated f c * and experimentally determined values, f e . As in the case 
of the mode-shape comparison, the computed coupled frequencies f c * generally agree 
well with the measured values f e for the first mode, but the results again compare less 
favorably at the higher modes. The deviation between f e and fc* for the first mode 
ranges from 0 to 15 percent, the greater number of these frequencies agreeing within 
8 percent. Agreement in the second mode frequencies was also within 8 percent except 
for a maximum deviation of 20 percent in one case. The third- mode coupled frequencies 
deviate from the measured values by as much as 20 percent except for a 35-percent dif- 
ference in measured and computed frequencies corresponding to the three mode shapes 
for the 5000 lbm (2268 kg), t = 0 yaw case already discussed. 

Longitudinal and tors ion modes, coupled analysis .- In addition to lateral response, 
the matrix Holzer analysis yields modes which exhibit dominantly longitudinal or torsion 
motion or solid- rocket- motor asymmetric yaw. These modes are shown qualitatively 
in figure 5 and are characterized as: 
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Mill 


(1) Solid-rocket-motor symmetric longitudinal motion combined with core longi- 
tudinal motion (fig. 5(c)) 

(2) Solid- rocket-motor asymmetric pitch motions combined with core torsion 
motion (fig. 5(d)) 

(3) Solid- rocket-motor asymmetric yaw motion combined with core longitudinal 
motion (fig. 5(e)) 

(4) Solid- rocket-motor asymmetric longitudinal motion combined with core yaw 
motion (fig. 5(f)). 

Experimental data were obtained only for the pitch-torsion modes. Symmetric longi- 
tudinal and asymmetric yaw motion was not observed during the tests since the experi- 
mental investigation was concerned primarily with lateral response of the three-body 
configuration. Consequently, comparison between measured and computed data is pos- 
sible only for the pitch-torsion mode. Nevertheless, for completeness, the longitudinal 
and yaw frequencies given by the matrix Holzer coupled mathematical model for the 
26 000 lbm (11 790 kg), t = 0 configuration are given in table IX. 

TABLE IX.- COMPUTED FREQUENCY FOR LONGITUDINAL, YAW, AND 
PITCH MODES OF SOLID ROCKET MOTOR 
(26 000 lbm (11 790 kg) payload] 


Mode 

Computed frequency, 
Hz, for - 

SRM asymmetric yaw; core longitudinally 

t = 0 sec 
13.7 

t = 105 sec 
44.3 

(fig. 5(e)) 



SRM asymmetric longitudinally; core 

17.3 

36.3 

yawing (fig. 5(f)) 



SRM symmetric longitudinally; core lon- 

23.9 

36.7 

gitudinally (fig. 5(c)) 



SRM asymmetric pitch; core in torsion 

15.3 

29.8 

(fig. 5(d)) 




The pitch-torsion mode was difficult to predict accurately for reasons somewhat 
similar to those given for the second pitch mode. A parametric study was made to deter- 
mine the effect of the aft outrigger torsional flexibility coefficient on the pitch-torsion 
mode. It was found that in order to obtain agreement with measured results, the aft out- 
rigger flexibility would have to approximate the flexibilities of the core and the forward 
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outriggers taken individually. By using recomputed coefficients, as indicated in table II, 
the results shown in figures 14 and 17 were obtained. 

Coupl ed analysis compared with uncoupled analysis .- The matrix Holzer analysis 
permits the assumption that combined pitch-torsion motion may be uncoupled from com- 
bined yaw -longitudinal motion. Vibration data were computed by using the uncoupled 
version of the analysis to determine whether the added complexity of coupling produces 
results in better agreement with the measured data. Uncoupled frequencies for repre- 
sentative configurations are given in table V. Negligible variations were noted in com- 
parable values of the coupled and uncoupled frequencies and indicated that the coupled 
and uncoupled versions of the analysis were equally successful in predicting the first 
three lateral modes of the model response. 

CONCLUDING REMARKS 

A matrix Holzer analysis employing a three-dimensional coupled mathematical 
model was used to compute the free vibration modes and frequencies of a three-body 
clustered launch vehicle. Computed data for the three lowest lateral pitch and yaw 
modes were obtained for propellant loading conditions of the model corresponding to 
flight times ranging from lift-off (time 0) to solid rocket burnout (time, 105 sec) and for 
payload masses corresponding to 5000 lbm (2268 kg), 26 000 lbm (11 790 kg), and 
45 000 lbm (20 410 kg) full scale. 

In general, the coupled mathematical model analysis adequately predicts the first 
three lateral vibration modes of the model configuration for all payload weights and pro- 
pellant loading conditions. Agreement between the computed and measured data is good 
for the fundamental mode but becomes less favorable at higher modes. 

The coupled system analytical data were computed by using aft outrigger attach- 
ment flexibilities reevaluated on the basis of experimental evidence. Computation of 
mode shapes and frequencies assuming a negligible restraint against rotation of the solid 
rocket motors in the pitch plane resulted in a decided improvement in agreement between 
the measured and computed pitch-mode response. This study indicates that correct 
treatment of intertank connections is critical for certain modes. 

Mode shapes and associated resonant frequencies computed by the coupled system 
analysis were compared with corresponding measured data and other computed data 
obtained from analyses using an uncoupled mathematical model. The results indicate that 
an uncoupled analysis is sufficient for this three-body configuration. 

Modes associated with dominantly longitudinal motion and yaw motion of the solid 
rocket motors were obtained analytically. Since no experimental data were available 
for comparison, the accuracy of these predictions was undetermined. 
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The use of dynamic models in an experimental test program in conjunction with an 
analytical study provides a valuable tool for accurate appraisal of the dynamic behavior 
of the physical vehicle. The experimental program provides an overall picture of the 
vehicle vibration characteristics and, in addition, can provide dynamic characteristics 
of the complex components of the model for which an accurate analysis is extremely dif- 
ficult or impossible. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., February 7, 1968, 

124-08-05-22-23. 
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APPENDIX 


TRANSFER MATRICES USED IN ANALYSIS 

The analytical procedure outlined in the main body of the present paper follows 
closely the analysis presented in section IC-1 of reference 6. The transfer matrices 
for the elemental massless rods and lumped masses employed herein differ somewhat 
from those presented in reference 6 and are therefore discussed in this appendix. Rela- 
tionships defining the outrigger flexibility coefficients are based on a sign convention 
which also differs from that of reference 6. Since these equations are used in the digital 
program written by the authors of reference 6, they are included in this appendix. 

Symbols 


a length of aft truss 

b eccentricity of aft truss, 2.6 in. (6.6 cm) 

EA(x) longitudinal -stiffness distribution 

EI(x) yaw-bending-stiffness distribution 

EI'(x) pitch-bending- stiffness distribution 

e eccentricity of mass center of gravity relative to cross-sectional centroid 

GJ(x) torsional- stiffness distribution 

I Q bending inertia of branch mass about its own center of gravity 

I^ bending inertia of lumped mass excluding flexible branch 

effective bending inertia of lumped mass including branch 
torsional inertia 

KAG(x) shear stiffness distribution 

K B rotational spring for branch mass 
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L 

l 

m Bi 

“Bi 

m i 

m li 

Pp?2 

r 

S B 

S 4>i 

Cl) 

W B 

Bn 

Matrices: 

O 

[] 

M 

H 

M 

M 


length of massless rod element 

distance from center line of core to center line of solid rocket motor (SRM) 

ith flexibly attached branch mass 

effective lumped mass including branch mass 

ith rigidly attached lumped mass 

mass effective for longitudinal motion 

aft outrigger joints for solid rocket motors 1 and 2, respectively 

overhang distance of branch mass 

first mas.s moment of branch mass 

first moment of ith mass excluding branch mass 

circular frequency of complete system 

rotational circular frequency of flexible branch mass 

rotational circular frequency of flexible branch mass with I Q = 0 

denotes column matrix 

denotes square or rectangular matrix 

flexibility matrix for pitch- torsion motion of forward outrigger 
flexibility matrix for yaw bending motion of forward outrigger 
flexibility matrix for pitch-torsion motion of aft outrigger 
flexibility matrix for yaw -longitudinal motion of core at aft outrigger 
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M 

M 

M 

O] 

0] 

W 

N 

M 

M 

Vectors: 


flexibility matrix for yaw- longitudinal motion of solid rocket motor and aft 
outrigger 

submatrix of elastic matrix for yaw bending 
submatrix of elastic matrix for pitch bending 
submatrix of elastic matrix for longitudinal deformation 
submatrix of elastic matrix for torsional deformation 
submatrix of mass matrix for yaw masses 
submatrix of mass matrix for pitch masses 
submatrix of mass matrix for longitudinal masses 
submatrix of mass matrix for torsional inertias 


— ► denotes force and linear displacement vectors 

— ►*- denotes moment and angular displacement vectors 

Transfer Matrix for Massless Rods 

The matrix which expresses the relationship between the internal loads and dis- 
placements at the ith and i + 1 stations across a massless rod is found from the static 
solution of the beam equations. The rods were assumed to have constant section proper- 
ties in reference 6. In the present paper, however, provisions were made to allow for 
variable cross-sectional properties between the lumped masses. The matrix elements 
as given are generalizations of those of equation (1) on page 10 of reference 6 and reduce 
to them for constant cross-sectional properties. The mass and elastic transfer matrices 
are related to the beam elements as shown in sketch 1. 
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By definition, 
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dx 

GJ(x) 


0 

1 


i 


1 



' L _dx_ 
0 AE(x) 


0 


X 


Transfer Matrices for Lumped Masses 

The elements of the mass transfer matrices in reference 6 were derived for con- 
figurations shown in sketch 2(a) where the ith mass mj may be comprised of one mass, 
rigidly attached to the main beam with a mass mg spring mounted to it. In this paper, 
the transfer matrices were modified to include configurations where mass mg might 
be flexibly attached to mi as indicated in sketch 2(b). The physical properties are 
obtained either theoretically, or as is often the case with complex mountings, from exper- 
imental data. Measured characteristics were employed for the 5000 lbm (2268 kg) pay- 
load analyses and were found to have a significant effect on both calculated frequencies 
and mode shapes. Expressions for the branch mass matrix elements are given. It 
should be noted that the branch need not be a lumped mass and spring but may be a multi- 
mass system representable with a single degree of freedom. It must nevertheless be 
attached to the main beam in a statically determinate manner. 



Sketch 2 
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Outrigger Flexibility Coefficients 

Transfer matrices representing the forward and aft outriggers (which are assumed 
to be massless) are obtained from flexibility relationships between loads and displace- 
ments. At the forward attachment which cannot transmit bending moments or longitudinal 
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Sketch 3 

forces, a static analysis employing the sign conventions indicated in the accompanying 
diagrams yields the following expressions: 


Y IV 1 " Y III 1 ■ 0 ID* 


S' 

-v’ 

IVl 

^IVi “ 

>-0]< 

■ TIV1 > 

T 

y iv 2 + y iii 2 ■ e m 2 


- y w. 

0 iv 2 ‘ e m 2 

V- J 


- t iv 2 


for pitch bending-torsion motion (see ref. 6, p. 22) and 



for yaw bending motion (see ref. 6, p. 23). Subscripts 1 and 2 denote solid rocket 
motors 1 and 2, respectively. 

At the aft outrigger, the following pitch-torsion load deflection relationships (eq. (14) 
in ref. 6, p. 28) are derived: 
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Yaw-longitudinal properties at the aft outrigger are based on the core flexibility and 
the truss- solid rocket motor flexibilities considered separately. This treatment is neces- 
sitated by the fact that zero yaw bending moment is transmitted at joints Pj and P 2 , 
whereas all other attachment joints are assumed to be fixed. From the following 
schematic, 


V,Y 


Sketch 5 
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consideration of core flexibility (see ref. 6, pp. 16-17 and pp. 30-31) leads to 


r Y~'' 






Y 

4 >' 

Y’ 

e 




OO 


VIII 
VII j 
VIIo 


and truss solid rocket motor flexibility (see eq. (15) in ref. 6, pp. 29-30) leads to 





It should be noted that transfer matrices for the forward and aft outriggers are 
never obtained explicitly. Certain relationships which would be contained therein are 
employed to obtain solutions for the intermediate unknowns of equation (4). 
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Figured- Variation of nondimensional |U(X)| with frequency parameter X. Curves determined from matrix Holzer analysis computed for 

coupled and uncoupled mathematical model. 
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Figure 12.- Measured and computed lateral pitch response using the modified aft outrigger flexibility data. 

t = 0; 26 000 lb (11 790 kg) payload. 
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Figure 13.- Measured and computed lateral yaw response using the modified aft outrigger flexibility data, 
t = 0; 26 000 lb (11 790 kg) payload. 
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Figure 14.- Measured and computed solid rocket motor asymmetric pitch response using the modified aft outrigger flexibility data. 

t = 0; 26 000 lb (11 790 kg) payload. 
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Figure 15.- Measured and computed lateral pitch response using the modified aft outrigger flexibility data, 
t = 105 seconds; 26 000 lb (11 790 kg) payload. 
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Figure 16.- Measured and computed lateral yaw response using the modified aft outrigger flexibility data, 
t = 105 seconds; 26 000 1b (11 790 kg) payload. 
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Figure 17.- Measured and computed solid rocket motor asymmetric pitch response using the modified aft outrigger flexibility data. 

t = 105 seconds; 26 000 lb (11 790 kg) payload. 
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